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Abstract 

We present the dark matter (DM) extension (by N) of the minimal supersymmetric standard 
model to give the recent trend of the high energy positron spectrum of the PAMEL A/HEAT 
experiments. If the trend survives by future experiments, the MSSM needs to be extended. Here, 
we minimally extend the MSSM with one more DM component N together with a heavy lepton E, 
and introduce the coupling eRE'^NR. This coupling naturally appears in the flipped SU(5) GUT 
models. The model contains the discrete symmetry Zq, and for some parameter ranges there result 
two DM components. For the MSSM fields, the conventional i?-parity, which is a subgroup of Zq, 
is preserved. We also present the needed parameter ranges of these additional particles. 
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The existence of dark matter (DM) at the 23% level of the closure density [l| is largely 
accepted by the observations of the flat rotation curves of the velocities of halo stars, the 
simulation of the bullet cluster collision, and gravitational lensing experiments. So, the 
identification of the cosmological DM is of the prime importance in particle physics and 
cosmology. If indeed the DM particles of O(IOO) GeV mass with the weak interaction 
strength are abundant in galaxies, high energy positrons, antiprotons and gamma rays from 

:he DM annihilation is confirmed. 



DM annihilation has been predicted for a long time . If t 
the very weakly interacting DM possibility [3| is ruled out 



The satellite PAMELA experiment has already started to probe an interesting SUSY 



parameter space. 
GeV up to 60 GeV 



"heir recent report on the high energy positron observation, above 10 
5|, has already spurred a great deal of attention G,!?]. In fact, the same 
trend has been noticed earlier in the balloon-borne HEAT experiment but with larger error 
bars Q]. The charactersistic of the PAMELA/HEAT data between 10-50 GeV is a slightly 
rising positron flux, e'^ /{e'^ + e~) ~ 0(0. 1). If confirmed by independent observations such 
as PEBS baloon experiments 9|] and the AMS-02 experiments lOj, the implication is tanta- 
mount in that the most popular LNx DM scenario of the minimal supersymmetric standard 
model (MSSM) maybe in jeopardy. Even though the possible astrophysical explanations 
have been presented in ll| , here we focus on the particle physics explanation. On the other 
. land, we note that the same PAMELA data does not show any significant anti-proton excess 

The lightest neutralino x (LNx) of the MSSM has been favored as the cold dark matter 
(CDM) candidate. If the CDM is composed of just one component Majorana fermion such 
as the neutralino LSP, due to the Fermi statistics two annihilating LSPs of Fig. [I](a), being 
the identical fermions, can come close if their spins are anti-parallel so that the angular 
momentum is zero, or in the s-wave state. Accordingly, when the neutralinos are annihilated 
to a fermion-antifermion (//) pair, the spins of / and / should be anti-parallel. For the 
helicity flipping, the amplitude should necessarily acquire a factor mj in the s-wave state. 
Thus, co-annihilation of neutralinos into c, 6, and t quark pairs, if it is open, would dominate 
over all other channels in the MSSM. In view of the conventional MSSM CDM scenario, 
therefore, the recent reported PAMELA data on the high energy positron excesses is quite 
embarrassing because of the above difficulty faced in the MSSM CDM scenario. 

This leads us to consider a minimal extension of the MSSM so as to keep its most desirable 
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property 'super symmetry'. In the framework of spin-| CDMs, therefore, we extend the 
MSSM minimally to include two CDM components. 

In the MSSM, if the LNx is bino denoted as Xi then the cosmologically favored bino 
density in the universe is possible in the co-annihilation region In this paper we treat 
the LNx as the bino just for the sake of a concrete discussion. If the LNx contains the bino 
as a fraction, the discussion should be modified accordingly. 

As commented above, two annihilating neutralinos of Fig. [Uj^a), being the identical 
fermions, can come close if their spins are anti-parallel so that the angular momentum 
is zero. If the final electron and positron, going out back to back, have the same helicity, 
then they can make up the angular momentum zero. One such possibility shown in Fig. 
[U^a) shows that it is highly suppressed because the bullet of Fig. [T]^a) carrying an SU{2)w 
quantum number fe{H^)m3/2 a high suppression factor. On the other hand, if the outgoing 
electron and positron carry the opposite helicities, their spin is one and by emitting a photon 
the three particles final state can make up angular momentum zero. However, in this case 
there is a coupling suppression of order aem/7r. This possibility of high energy positron plus 
photon has been suggested in Ref. where a large enhancement factor of order 10^ is 
needed. A scalar DM such as sneutrino LNx will have the same fate as the bino LNx in this 
regard. 




(a) (b) 

FIG. 1: The bino-like neutralino annihilation, (a) The bullet carries an SU{2)w quantum number, 
(h) Here, the bullet can carry angular momentum 1. The helicities of electrons and positrons are 
shown by thick arrow lines. 

One may consider a very heavy neutralino dominantly producing W^W~ so that sub- 
sequent W decays provide the positrons. This requires a very large LNx mass 0(10 TeV) 
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which belongs to the neutral component of weak isospin 7 = 5 multiplet to satisfy sev- 



eral constraints as a DM candidate 13|. This is also not very minimal in the sense that 
many component fields are introduced and SUSY particle masses are raised substantially. 
In addition, the intermediate W^W~ will give some excess anti-protons also. 

This consideration presents a useful direction for constructing DM models with the ris- 
ing high energy positron spectrum above 10 GeV. Obviously, a Dirac particle DM can be 
considered Q Also a sp,„-l DM fl ca„ overcome the bino-like LNx difficulty, but here 
the spin-1 DM [7|] would give some excess anti-protons also. 

Thus, we extend the MSSM so that the DM annihilation produces high energy positrons 
but not excess anti-protons, thus providing a possible explanation of the PAMELA/HEAT 
data. This needs a special treatment of color singlet particles compared to quarks. So 
let us consider the next MSSM 



151 ] extended by a neutral singlet (A^dmMSSM). A^ is 
split into two Weyl spinors {Nji.Ni}. We also introduce an SU{2) singlet charged lepton 
E (= {Eji, El}). Without introducing this singlet lepton, we cannot achieve the goal of 
producing excess positrons: This charged lepton makes A^ stable in some region of the 
parameter space as we will discuss below. 

This A^dmMSSM seems to be very simple in the sense that new particles needed beyond 
the MSSM are minimal, just A^ and E. To explain the correct order of DM density and 
the production cross section of positrons, the weak interaction nature of WIMP must be 
kept. Let us introduce the continuous /^-symmetry so that the additional DM component 
A^ remains light down to low energies. The weak hypercharge Y and R charge of the singlet 
fields are 

Superfields : cr Nr N% Er E% 

Y : -10 0-1+1 (1) 

R : +1 i -f -| +1 

The R charges of the MSSM fields are as usual: the quark and lepton superfields carry 1 
and the Higgs superfields carry 0. The R symmetry allows the superpotential, 

W = /crE^Nr + hNl (2) 

where / and h are coupling constants, but the mass terms of A^ and E cannot be present 



in W . However, via the Giudice-Masiero mechanism [l6(], supergravity effects can generate 
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the fermion masses if the F term of a singlet 5* is developed, 

I d'e^ {XEnE% + >!NnN% + h.c.) . (3) 

Similarly the MSSM //-term is also generated. These masses are assumed to be of order 
the gravitino mass 777,3/2- does not develop a VEV. As usual, with the separate lepton 
number conservation the process yU — > 67 is forbidden. Without the separate lepton number 
conservation, the f/(l)_R allows i?^ couplings to and tr by couplings /' and /", in which 
case we need |/'| < 10^^. The /" bound is much weaker. 

Not introducing couplings to quarks, we will not introduce excess anti-protons. From 



string compactification with the doublet-triplet splitting 



for this purpose [18 1. 



ITj, the flipped SU (5) is best suited 



The mass terms induced by SUSY breaking can be written as an effective superpotential, 

W D my2ERE% + tti'./.NrN^^, (4) 

and soft SUSY breaking A and B terms are 

C 3 m3/2[afeRE%Nn + ahNl] 

+ml^[hERE% + + h.c, (5) 

where a and h denote dimensionless couplings. Eq. ([5]) violates the continuous R symmetry, 
leading to a discrete Zg symmetry. The i?-parity (Z2) is a part of this Zq symmetry. For 
the fields in Eq. ([T]), the (i?— parity, Zq) charges are 

eR(-,3), Nr{+,2), iV^(+,4), ^^(-,5), ^^(-,1), 

e/?(+,0), iVH(-,5), Ar]^(-,1), Er{+,2), ^^(+,4). (6) 

The SM fields are given Zq = and their superpartners are given Zq = 3, which is just the 
i?-parity for this MSSM subset fields in our A'dmMSSM. The proton longevity is protected 
by the U{1)r. The U{1)r forbids URdRdR and the dimension 5 operators qR^RgRlR and 
URURdRCR which are allowed by the conventional -R-parity alone. If SUSY breaking induce 
them in the A^dmMSSM, they must be highly suppressed since there are no simple diagrams 
for them. Because of the -R-parity of A^, A^ cannot be a candidate for the singlet heavy 
neutrino of the seesaw mechanism. 
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We intend to introduce two stable particles, one the LNx and the other the lightest 
matter particle (LMP). For the superfields N^i + NfiO and Ef> + Er9, we assume > rriN, 
> ruE, and rriE > niiq. We assume that is lighter than cr and Er, and hence TV 
is taken as the LMP. The LMP carries = 5 which cannot be made by SM particles 
(carrying = 0) alone. If x is much heavier than A^, then the decay, x ~^ SA^'^ + + e~ 
suppressed by e^, E'^, and Nr propagators, is possible.^ Below, we consider the case of A^ 
being lighter than x- The decay rate is estimated as, F ~ g''^\f\^\h\'^ / M^Mjm\. But 
we require < Sm^r for a successful two DM components in the universe. In this way, we 
have two DM components, x ^'^^ ^ ■ Now, we take the bino as the LN^ and the fermionic 
partner A^ as the LMP. 



FIG. 2: A typical diagram for a bibo x annihilation. NuNf^ annihilation to e^e is also possible. 

For this idea of two DM components to work, we must satisfy the following: 

• The annihilation through Fig. [2] should be allowed. Namely, + rriN > uie + rrie. 

• X OT N should not decay by the diagram Fig. O and < niN + niE + rrte should be 
satisfied for > rriN. The case niisf > turns out to be impossible. 

• The E decay is allowed by the interaction E^ ^ cr+Nr, requiring niE > Mff+nie. 

• The bino decay x ~^ SA^'^ + + e~ is forbidden kinematically, < 3mj\f + 2me. 

^ Consider the Feynman diagram, x ^ Ee^N'^ihy e propagator) (e^iV)(by E decay) + e^N^' + 



en 



X 



} 




6 




rriN 



rriN 



rriE 



FIG. 3: In the — tue plane, the kinematically allowed mass region is shaded for a typical 
mass value of rriN- For Mjy > 2mN, N — > NN decay is possible. For > SniN, the decay 
X — > 3N^e~^e~ is possible. 



By the interaction (Ej), E'^ by Nr propagator {E'^ ^l^r ~^ <^l^R^r) directly 
(A?"^ 2Nr by ([2])) decay to 2A''R + e^ and 2Nr, respectively. For the latter two body decay 
the total decay rate is 



R) 



IGvr 



2m|r 

N 



(7) 



where /3 



N 



[I 



m 



Nl 



/4M^)^/^. The mass limit niE > 102.5 GeV for a stable heavy lepton 



E [19| does not apply here. But for our unstable E, let us use this bound as a guideline. 
For the LNx mass we take the usual MSSM estimate of order 100 GeV. Here, we assume 
that masses of and Nr^l are relatively small such that the needed kinematics are satisfied. 
Neglecting the electron mass, in Fig. [3] we plot the allowed region in the —mE plane for 
a specific rriN. As far as the decay rate is large enough so that E and A^ decayed before 1 
second, these decays are not problematic in the nucleosynthesis. Note that given an allowed 
phase space this is easily satisfied with not too small couplings / and h. 

If X and A^ each constitutes 50% of the CDM density, the annihilation diagrams of Fig. 



[2] account for | of possible encounters of x ^^id A^: xXi NN, xA^, and Nx- Among these, the 
dominant contributions come from x^R ^R^r annihilation. In principle, the ratio of 
X and A^ abundance is determined if high energy dynamics is completely known. 

Using the interaction ([2]), the cross section da{xN — > E~)/dQ of Fig. [21 in the center 
of momentum frame with the incident three momentum p, is calculated in the small |p| and 
the large limit as 

da |^7P 1 m^mN^/s[3%{l + niE/ ^/sY 



dn 



1287r2 IpI [^2 _ ^2 + m^y^p%\ 
2^Pl 1 



xU + 



PI 



p|cos^k (8) 



where Q is the angle between three momenta of x and e+, and = 1 — m\/ s with ^/s = 
+ m]\f. [We included the A^A^'^ annihilation process also, which however is suppressed for 
a large Eji mass.] A similar expression holds for the charge conjugated final states. Then, 
the velocity averaged cross section is calculated as 

In Fig. m we present our estimate of the positron excess for typical values of = 200 
GeV, niN = 80 GeV, = 200 GeV, = 400 GeV, and Mg = 220 GeV (thick green 
line), 250 GeV (blue dash line), 280 GeV (brown dash line). For a good fit, we need a small 
difference for Mg — M^. 

To compare with the observations, basically we use the astrophysical background flux 
given by ^fi^ = 4.5^°-V(l + 650^2.3 ^ isooE^^ and ^fj = 0.1QE~^-^/{l + UE^-^ + 



3.2^2-15) + 0.70^0-7(1 + 110^1-5 + 580E^-2) 
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221. The deviation of the PAMELA data 



rom this curve at low energy(< 10 GeV) can be explained by the solar modulation effect 
The calculation of the positron flux from a given particle physics model is well described 
261 . The positron flux is given by $e+ = Ve+^/'in, where Vf.+ is the velocity 



in Refs. 
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of the positron and ^ is the positron number density per unit energy, ^ = dN^+fdE. ^ is 
determined by the diffusion-loss equation using the various cosmic ray data as described in 



2l|. Under the steady state approximation, the solution of the diffusion-loss equation is 
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FIG. 4: The positron fraction from our model with My = 200 GeV, = 80 GeV, rriE 



200 



GeV, = 400 GeV and Mg = 220 GeV (thick green line) and B = 7. Mg for 250 GeV (blue 
dash line) and 280 GeV (brown dash line) are also shown. The pink band is the positron fraction 
coming from and NN annihilations and the green band is this positron excess on top of the 



astrophysical background (the thick darkblue dash line) ]20 . \2i l. The width of the band shows the 
uncertainty from the positron propagation model. The PAMELA data are the red dots Q/, and the 



various small dots represent the observed positron cosmic ray data 



given by a semi-exact form 



where I{\d) is the halo function which has the halo model dependence but is independent 
from particle physics and i? > 1 is a possible boost factor coming from the DM halo 
substructure |21 |. 

In conclusion, it is likely that the LNx of the MSSM cannot explain the high energy 
positron spectrum of the PAMELA/HEAT data [sl, O], simply because of the angular mo- 
mentum constraint. So, we extended the MSSM to keep low energy supersymmetry with 
additional fields and E with the U {1)r symmetry, to obtain two DM particles x and A^. It 



is the minimal extension toward explaining the PAMELA/HEAT data. In the A^dmMSSM, 
it is shown that a wide range of the parameter space for two DM components is allowed. 
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